function [k_c, Ma_c, alpha0] = growth_compare

p = params;

% initial guess for k_c is 1e-1
k0 = 1e-1;

% should start with Ma = 1e6;
[k_c, Ma_c] = weighted_neutral_stab_curves(k0);

alpha0 = zeros(size(k_c));

fprintf('-------------------------\n');
fprintf('computing alpha...\n');
fprintf('k_c\t\tMa_c\t\talpha_0\n');
fprintf('-----------------------------------------------\n');
for n = 1:length(k_c)
    p.Ma = Ma_c(n);
    [t, tmp] = time_eigs(k_c(n), p);
    alpha0(n) = tmp(1);
    fprintf('%.3e\t%.3e\t%.3e\n',k_c(n), Ma_c(n), alpha0(n));
end